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ABSTRACT 

The prime evidence underpinning the standard ACDM cosmological model is the 
CMB power spectrum as observed by WMAP and other microwave experiments. But 
Sawangwit & Shanks have recently shown that the WMAP CMB power spectrum 
is highly sensitive to the beam profile of the WMAP telescope. Here, we use the 
source catalogue from the Planck Early Data Release to test further the WMAP beam 
profiles. We confirm that stacked beam profiles at Q, V and particularly at W appear 
wider than expected when compared to the Jupiter beam normalised either directly to 
the radio source profiles or using Planck fiuxes. The same result is also found based on 
NVSS and WMAP 'CMB-free' source catalogues. Further, the WMAP source fluxes 
demonstrate a non-linear relation with Planck fluxes, as previously found between 
WMAP and ground-based fluxes. Importantly, we find that applying this non- linear 
relation to the Jupiter beam profile results in an excellent fit to the observed radio 
source profiles. Also, stacked SZ decrements of ~ 160 galaxy clusters observed by 
Planck appear to be underestimated in the WMAP data by a factor of Ri 2 in the 
Q, V and W bands. The strength of this discrepancy at Q means that it cannot be 
fully explained by the beam profile problem and thus the WMAP SZ discrepancy 
remains unexplained. Beam profile systematics can have significant effects on both 
the amplitude and position of the first acoustic peak. In particular, a wider beam can 
move the position of the first peak to significantly larger wavenumbers with potentially 
important implications for cosmology. 
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1 INTRODUCTION 

CMB experiments such as WMAP have made significant 
progress in the study of the primordial temperature fluctu- 
ations. Their best fitting power spectra strongly support a 
spatially flat, ACDM, universe. This model requires rela- 
tively few parameters, yet apparently manages a compelling 
concordance between a variety of other cosmological data; 
SNIa, Large Scale Structure and Big Bang Nucleosynthe- 
sis. Although the statistical errors on these power spectra 
are small, this precision does not necessarily imply accuracy 
and there remains the potential for systematic errors to alter 
these conclusions. 

Indeed, several anomalies between ACDM and the 
WMAP data have been discussed. Typically these 
have involved the l arge-scale temperature multipoles eg: 
IjBennett et al.ll201ll ). However, other anomalies in the CMB 
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at smaller scales have also b een detected, connected in p ar 



ticular with radio sources (|Sawangwit fc Shankd 2010al 
and SZ dec rements from galaxy clusters (JMyers et al.ll2004l : 
iBielbv fc Sh anks 2007) 

Radio sources are sometimes regarded as a contam- 
inant in CMB temperature maps. However, radio point 
sources prove particularly interesting because they pro- 
vide a complementary check of the beam measured by the 
WMAP team fr om observations of Jupiter (IPage et al.ll2003l : 
iHillet al.|[2009l '). Jupiter has a flux of ^ 1200 Jy which is 
« 3 orders of magnitude higher than radio source fluxes or 
CMB fluctuations. This high flux has advantages in terms 
of defining the wings of the beam profile but is a disadvan- 
tage in that the calibrating source is much brighter than 
typical CMB fluctuations. Furthermore, Jupiter only checks 
the beam o n the ecliptic whereas radio sources are spread 
over the skv. lSawangwit fc Shankd (|2010a bl made a stacked 
analysis of radio point sources and found evidence for a 
wider beam than WMAP measured using Jupiter. A tenta- 
tive detection of a non-linear relation between WMAP fluxes 
and ground based radio telescope fluxes was also found. A 
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thorough analysis of possible systematics did not find an ex- 
planation and we return to these issues later in this paper. 
The beam profile of a CMB telescope like WMAP is crit- 
ical because it smoothes the temperature anisotropies and 
therefore needs to be known accurately to produce the final 
power spectrum from temperature maps (|Page et alj 120031 : 
iHill et al.llioogi '). 

Various authors have noted small-scale anomalies with 
respect to the SZ decrements measured by WMAP. SZ 
decrements are created when the CMB inver se Compton 
scatte rs off hot electrons in galaxy clusters. iMvers et al.l 
((2OOJ) first stacked WMAP data at the positions of galaxy 
clusters and suggested that the profiles wer e mor e ex- 
tended than expected. iLieu. Mittaz. fc Zhang (|2006l ) and 
iBielbv fc Shanks! (|2007l ') then found that the SZ decrements 
from WMAP were reduced compared to X-ray predictions, 
possibl y due to the WMAP b eam being wider than ex- 
pected. HieMAlEanki (|2Q03) also found that the WMAP 
decrements were 3-5a lower tha n the ground-based SZ mea- 
surements bv lBonamente et al.l (|2006) in 38 X-ray luminous 
clusters. These anomalous WMAP SZ res ults have been con- 
firme d by a variety of re-analyses, (e.g. iDiego fc Partridge! 
!2009l). including a st acking analysis by the WMAP team 
iKomatsu et al.l!201l] '). 

In their ESZ sample, the Planck team find excel- 
lent agreement with the se lf-similar X-ray estima tes of the 
SZ decrement (jPlanck Collaboration et al.l !201ll ). This is 
corroborated by the ground based South Pole Telescope 
Collaboration with their blind SZ selected cluster sample 
iMroczkowski et al! !2009! 'l. This co mpounds the question 
of wh y W MAP SZ analyses frorn !Lieu. Mittaz. fc Zhanel 
l|2006l ) and !Bielbv fc Shanks! J2007!) fa i led to find such an 
agreement. Although !Komatsu et al.l (|201l! ') claimed the 
WMAP discrepancy was only a 1.5-2a effect, they focussed 
on comparing the integrated Y parameter rather than the 
inner SZ profile which generally shows a more significant 
discrepancy (Bielbv fc Shank&.2007l '). 

In this paper we use the recent Planck Early Data Re- 
lease and other radio source data to re-investigate both the 
WMAP radio source beam profile and SZ anomalies. The 
Planck Early Release Compact Source catalogue (ERCSC) 
is of particular interest and provides the basic parameters 
of radio sources and SZ clusters from the Planck CMB 
maps. Although, the corresponding temperature maps from 
which these were estimated have not been released, both 
radio source fiuxes and SZ profile parameters are available 
as measured by Planck. We can therefore use these to com- 
pare WMAP and Planck radio source fluxes directly and 
also to make WMAP stacks centred now on the new radio 
source and SZ cluster lists from Planck. From these stacks, 
the WMAP beam proflle can be inferred and the SZ results 
from WMAP and Planck compared. Given the higher angu- 
lar resolution, lower noise and different calibration strategy 
for Planck, this comparison will allow new insight into the 
robustness of the WMAP CMB analysis. 



2 DATA 

2.1 Planck Early Data Release 

The Planck team have recently made their first release of 
data collected by the Planck satellite between 13 August 



2009 and 6 June 2010 (amounting to ~ 1.5 full sky sur- 
veys). This early data release is concerned solely with the 
foreground contamination in the CMB maps. The two sets 
of catalogues relevant to this paper form the Early Release 
Compact Source catalogue (ERCSC). These are the Radio 
Source catalogues and the SZ catalogue. 

2.1.1 Planck Radio Sources 

The Early Release Compact Source Catalogue (ERCSC) 
lists all the high reliability radio sources with accurate flux 
determinations. The ERCSC has been quality controlled so 
that ^ 90% of the reported sources are reliable, > 5a detec- 
tions and that the fluxes are determined within ^30% accu- 
racy. The catalogues are band specific and for the bands of 
interest (1/ ^ lOOGHz) are created using the 'PowellSnakes' 
method, a Bayesian multi-frequency algorithm for detecting 
discrete objects in a random background. Flux estimates 
were obtained by use of aperture photometry within a cir- 
cle of the beam's FWHM. For the case of unresolved and 
potentially faint point sources, the Planck team recommend 
the use of the parameter FLUX and its corresp onding error, 
FLUX_ERR iPlanck Collaboration erall!201l! ). 

We reject any extended objects from the catalogue 
to maintain an unresolved sample with which to test the 
WMAP beam profile. To do this we make use of the Planck 
quality tag 'EXTENDED'. This is defined by comparing the 
source areal profile with the 2-D Planck beam. An addi- 
tional quality fiag 'CMBSUBTRACT' has also been pro- 
vided, which reflects on the quality of the source detection 
in a map with the best estimate of the CMB removed. In this 
paper we reject an y sources with a value of 2, as the Planck 
team recommend, (!Planck Collaboration et al.l!201ll ). 

We have also rejected sources within 4° of the LMC and 
sources at low galactic latitude, |6| < 5°. We also rejected 
sources flagged with high astrometric error and any sources 
with WMAP counterparts that are flagged as contaminated. 

Band corrections between WMAP and Planck have 
been ignored in the absence of information on the 
spectral index. This factor i s in any case small 
jPlanck Collaboration et al.1 !20lil '). The fuU details of the 
catalogue construction and comp osition are described by 
\Planck Collaboration et al.l (|201l! ) and briefy overviewed in 
Table [2^ The number of sources, A'^, selected in each band- 
pass are also shown in this Table. 

2.1.2 Planck SZ Catalogue 

The Early SZ (ESZ) catalogue lists all the robust and ex- 
tensively verified SZ detections in the first data r elease . 
As described by !Melin. Bartlett. fc Delabrouilld l!2006l ). 
the Planck team extract the integrated Y parameter 
using a Multifrequency Matched filter (MMF3) method 
{Planck Collaboration et al.l |201l|). The algorithm is run 
blindly on all-sky maps, assuming the characteristic SZ spec- 
tral signature and self-similar cluster profile. 

In the Early Release of the Planck SZ catalogue, only 
data from the lOOGHz frequency channel or higher has been 
used to study the SZ effect. This is to avoid the deterimental 
effect on S/N from beam dilution caused by the larger beam 
sizes of the lower frequency channels. At the higher frequen- 
cies, the Planck beam FWHM is typically ^ 4.'5. The fuU 
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Freq FWHM A^ Flux Limit 
[GHz] (') [Jy] 



100 
70 

44 



9.94 
13.01 
27.00 



251 
134 
114 



0.344 
0.481 
0.781 



Table 1. Summary of the Planck bandpass parameters and the 
flux range of the sample we use from the ERCSC after w e have 
taken our quality cuts, l lPlanck Collaboration et al.ll201 j) . N is 
the number of radio sources to the specified limit. 



details of t he catalogue construction and com position are 
described m \Planck Collaboration et al.l l|201ll '). 

The catalogue provides estimates of the SZ flux, extent, 
redshift and position. It consists of 189 clusters, all detected 
at high S/N {^ 6) with 95% reliability. Whilst the sample 
is primarily composed of known clusters (169/189), it pro- 
vides a wealth of new information as it gives the first SZ 
measurements for ~ 80% of the clusters. In this paper we 
only make use of clusters with known redshifts and these all 
have X-ray observations. We therefore consider 159 clusters, 
including Coma. For this sample the redshift range spans 
z G [0.0126, 0.546] with a mean redshift of 2 = 0.179. 

2.2 WMAP Data 

We will be using the 7-year WMAP temperature maps ob- 
tained from the LAMBDA CMB resource. We work with 
the Nside ~ 512 Healpix maps resulting in a pixel scale of 
7'. We use the foreground unsubtracted temperature band 
maps for Q,V and W. We estimate the Jupiter beam in each 
band by averaging the 7-year beam profiles from the various 
DA's. To avoid Galactic synchrotron contamination we have 
used the extended temperature mask which admits 71% of 
the sky. When working with radio point sources we instead 
use the point source catalogue mask. 

We have used the 7-year WMAP point source cata- 
logue. These sources are detected at least the 5a level in 
one WMAP band. For a flux density to be stated, the 
detecti on must be above th e 2a lev el in that band. Fol- 
lowing ISawangwit fc Shanksl l|2010ah we ensure that the 
sources are genuinely point source by matching to the 5GHz 
(~ 4'. 6 resolution) catal ogues from the Greenbank North- 
ern sky Survey fGB6). ' Gregory et al.l lll996r) . or Parkes- 
MIT-NRAO (PMN'l. iGrrffith fc Wrightl(| 19931 ), surveys. The 
WMAP te am also provide a 7 yr CMB-free catalogue as de- 
scribed bv lCold et al.l (|201ll ). This catalogue has been cre- 
ated with the objective of detecting point sources free of 
boosting by CMB fluctuations. After the 5GHz matching, 
we proceed with a raw catalogue with 471 sources and a 
CMB-free catalogue with 417 sources. 



3 PLANCK RADIO SOURCE FLUXES AND SZ 
CLUSTER DECREMENTS 

3.1 Conversion of Radio Flux to Temperature 
Profiles 

The Planck ERCSC provides us with the source flux den- 
sity, error and a few parameters on the source characterstics 



Band 


Freq 
[GHz] 


FWHM 
(') 


(sr) 


r/7 

ll^KJy-^ 


9H 


W 
V 

Q 


94 
61 
41 


12.6 
19.8 
29.4 


2.097-10-5 
4.202-10-5 
8.978-10-5 


179.3 
208.6 
216.6 


1.245 
1.099 
1.044 



Tabl e 2. Summ ary of th e WMAP bandpass parameters taken 
from Hill et al.l (2009i ) and ljarosik et al.l l|201ll ). See text for def- 



and detection. To enable us to translate the Planck fluxes 
into WMAP observables we need to convert the source flux 
density into an observed peak Rayleigh- Jeans t emperature 
using the conversion factor, F (|Page et al.ll2003l ). 

ATflj(O) = StotT (1) 

A simplified form for F can be found. IPage et al.l (120031 ). 



-^/7 



ZiK,i)Vf, i Lbean 



(2) 



Here Ve is the effective frequency and F-^-^ denotes the fact 
that the majority of the WMAP sources have a power index 
a ~ —0.1, approximately that of free- free emission. 

The WMAP temperature maps are given in terms of 
the thermodynamic temperature. At the WMAP frequencies 
and CMB temperature, the Rayleigh-Jeans temperature is 
appreciably different from this. We therefore correct between 
the two temperature differences, using eq©, where x' — 
hi/fktTcmb and Tcmt = 2.725K is the monopole temperature 
of the CMB. 



ATt 



g{iy)ATRj. 



ATr 



(3) 



The observed WMAP temperature profiles therefore take 
the form. 



AT{9) = AT(0)&"(6') 

^g{v)vffStoth%e). 



(4) 



We see the beam dependence of the observed profile is 
twofold. The shape is dependent on the symmetrized beam 
profile b'iO) (normalised to unity at 6* = 0°), while the scale 
is normalised by the beam solid angle associated with F-'-' . 
A summary of the assumed values of g{v) and F-^-^ are pro- 
vided in Table [2] 



3.2 Planck SZ Decrements 

Planck presents its observed decrements using an SZ model 
fit parameterised by the total SZ signal within the cluster ex- 
tent. Here we briefly describe this model so that the Planck 
results can be compared to the stacked WMAP temperature 
decrements. 

Clusters are significant reservoirs of gas which will result 
in a SZ distortion to the CMB described by the Compton y 
parameter, 
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AT{e) = n^tj{x')y{e) 



(5) 



j(x') is the spectral function, where x' = hi/fkbTcmb (see e.g. 
iRefregier. Spergel. fc Herbig l200d ) The integrated Y pa- 
rameter is the total SZ signal, which is simply the integration 
of the Compton y parameter on the sky. 



Type 


Po 


C500 


7 


a 


/3 


AllTittcd 


8.403(5^)-i 


1.177 


0.3081 


1.0510 


5.4905 



Table 3. Summary of the Planck NFW pa rameters as used in 
eqn. llOl and described bv lArnaud et al.l 1120101) . These are the same 
parameters as used by the Planck team, the AlhFitted set. 



F= f ydQ. (6) 

Alternatively, if we integrate over the cluster volume. 



F = 



o-t 



/ 



PdV. 



(7) 



However, we are observing a 2-D projection of the clusteip 
on the sky. The angle 6 we observe on the sky, corresponds 
in 3-D to a cylindrical bore of the cluster of radius R = 
9Da{z). where Da is the angular diameter distance. The 
observed integrated Y parameter therefore takes the form 
IJArnaud et al.ll2010l ). 



Ycyl{R) - 



R pRtot 2P{r')r'dr' 

2-Krdr ' 



1^12 _ j.2')l/2 



(8) 



— Yaph{Rtot) 



Rto 



(9) 

To predict the SZ effect implied by eq((9| we have to 
make a choice of the pressure profile P{r). Historically it 
has been common to fit the SZ profile with an isothermal /? 
model, ([Cavalier e fc Fusco-F omiano 1976). However, X-ray 
observations have shown that the assumption o f an isother- 
mal g as breaks down at th e cluster outskirts, (jPratt et al.l 
I2OO7I : IPiffaretti et all 120051). To account for this additional 
complexity, iNagai, Vikhlinin. fc Kravtsovl l|2007 'l therefore 
proposed using a Generalised NFW (GNFW) profile for the 
pressure instead. The profile is scale invariant in that it is 
independent of absolute distances and is instead a function 
of the dimensionless scale x = R/Rsoo- The profile takes the 
form. 



V{x) = 



Po 



(C5003;)T[1 + (csooa;)"']"'-^'''" 



(10) 



where V{x) — P M/P^qq where P500 i s the characteristic 
pressure defined bv lArnaud et al.l ([20101). 

Here we have a five parameter fit to the pressure pro- 
filqj [^0, C5oo,7, ct, /^l- This allows independent specification 
of the pressure in the cluster core (7), main-body (a) and 
outskirts (/3). In Table [3] we outline th e para meters used 
by Planck, as taken from lArnaud et al.l ((20101) . The char- 
acteristic parameters of the cluster are M500, P500, Rsoo 
(see Appendix [A]) where the 500 denotes the fact they 
are evaluated within the region where the mean mass den- 
sity is 500 times greater that the critical density pcrit{z). 



^ The cluster is assumed to be spherical. 

^ Note the (B introduced here is different from the /3 associated 

with the isothermal fit to the pressure profile. 



The Planck team extract the integrated Y parameter us- 
ing the Multifrequency Matched Filter (MMF3) method 
(jPtoracfc Collaboration et al.ll201ll ) based on the above self- 
similar model. The integration is done to the angular clus- 
ter extent corresponding to 5-R500, which they also report 
(^5iJ50o)- Their errors on the SZ integrated Y and angular 
extent, ^5^500, are estimated by iterative re-measurements. 

3.3 SZ Temperature Profile Reconstruction 

We now proceed to invert the Planck data to provide us 
with expected WMAP temperature profiles. (See Appendix 
1X1 for the details of this derivation) . /.From the Planck val- 
ues for y(5i?5oo jj and Sshsoo, and using J{x) and I{x), the 
cylindrical and spherical SZ templates (see egns. ET)) we first 
obtain Ycyi{R) via eqn. lAQI 



Y,yi{R) = Y,yi{5R50o)(l - J 



J{x)_ 
(5) 



(11) 



From this integrated Ycyi{R — 9 ■ Da(z)), we want to derive 
the angular dependence of the Compton y parameter, where 
yiO) = ^Y^yiiO), and so 



//IN Ycyi{5R5oo) d , . ,s 

y^'^ = - 1(5) rfn(-^("))- 



(12) 



The self-similar model therefore predicts an SZ temperature 
decrement. 



^Tsz{e) = nmtjix') 



Ycyl{5R5' 



1(5) 



dn 



{J{x)) 



(13) 



where Ycyi{5R5oo) is the integrated Y given in ESZ. 

We have corroborated this method with an alternative 
calculation of Ycyi{R) which is independent of the Planck 
provided y(5i?5oo) and instead solely uses Osrsoo- The strat- 
egy is to eva l uate M500, to find P500 using equn. 13 of 
lArnaud et al.l l|2010f) and hence Ysoo- Ycyi{R) is the prod- 
uct of this and the cyclindrical scaling function J{x) (see 
ean. lA5|l . We find that the inferred SZ profiles between our 
two methods agree very well with < 10% difference at all 
scales with this small difference being caused by the sec- 
ond method's dependence on the Y — M calibration and the 
specification of the cluster physics. 



3.4 Isotliermal profiles from SZ parameters 

We have seen that the reconstructed Planck decrements are 
actually reasonably insensitive to the SZ template, at least 
within Psoo- This even applies when the self-similar model 



It is assumed that YcyiibRsoo) = ^(S-Rsoo) 
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is re placed by an i s other mal model. On the basis of Fig. 
11 of lArnaud et al.l ((20101) which compares Y values as es- 
timated from the REXCESS sample, the isothermal model 
requires a "^5^500 parameter « 3 x larger than the self-similar 
model but the ^^500 values are more similar. This gives us 
a route to produce an SZ decrement from the isothermal 
model based on the measured Planck Y value. The isother- 
mal model then has a Compton-y parameter of the form, 



y(&) = yo 



1 + 



(14) 



Equating Yp{5R5oo) = 3.2 x YsR^gg, we solve for the ampli- 
tude yo using Yfi{5RBoo) = Jj^"*^-^™"^ y{e)dn. We find that 



yo = Y,3 



2 ''' 2 



1 + 



-M- 



(15) 



We can therefore predict the SZ t emperature decrement 
using eq([5| and eq(ll4l). Fol l owing iLieu. Mittaz. fc Zhang 
l|2006t ) and iBielbv fc Shanks! l|2007h we assume 9c = 9.8' 
and P — 0.75 from the Coma cluster as typical isothermal 
model parameters, scaling Oc with redshift using a f^A ~ 0.7, 
fim = 0.3 angular distance relation. 

Clearly, the isothermal model may represent the Planck 
results less accurately than the self-similar model from eqn 
1131 But we shall see that, for fixed integral Y, the isother- 
mal model produces lower central temperatures and it will 
be interesting to see if that helps the agreement with the 
WMAP stacked SZ decrements. 



4 CROSS-CORRELATION METHODS 

Our cross-correlation/stacking procedures for measuring 
both radio poi nt source profiles and SZ decrements a r e sim - 
ilar to those of lMvers et all ll2004l).[Bielbv fc Shanksl (|2007l ') 
and then as updated bv lSawangwit fc Shanksl (2010a|). Ulti- 
mately, we shall be stacking/cross-correlating WMAP data 
around radio source positions and cluster centres from cat- 
alogues, particularly from Planck ERCSC. To estimate a 
decrement for an individual source we use, 



W(a) n PLANCK'WMAP (82) 
ATCA'WMAP (6) 
" IRAM*WMAP (43) 



Q 10^ 





1'* slope = 0.70±0.04 

T i, inlercepl = 0.13±0.03 



1 10 

PLANCK S,„„a„ / ATCA S,-^,^ / IRAM S,„c;h. [JyI 



V(b) n PLANCK'WMAP (83) 



£ 10 ^ 



S B 





slope = 0.83±0.04 
intercept = 0.07±0.04 



1 10 

PLANCK S,„„„[Jy] 



- Q(c) n PLANCK'WMAP (76) 
Q 10^ 




1 10 

PLANCK S^^;,^[Jy] 



Figure 1. A comparison between the WMAP7 fluxes, Planck 
and ground based source fluxes. 



experimented with both alternative sub-fields and methods 
such as bootstrap resampling, finding approximately equiv- 
alent results. 



ATi9) = Yl 



T,{e)-T 

m{9) 



(16) 



where the sum is over the pixels, denoted i, within a cir- 
cular annulus centered on the angular separation 9. Ui rep- 
resents the number of pixels within the annulus and Ti is 
the temperature recorded for the pixel i. T is the average 
background temperature which can either be estimated lo- 
cally in a sur rounding annulus in a 'photo metric method' or 
globally fsee ISawangwit fc Shankal2010al 'l. These two back- 
ground estimates make no difference in the stacked results 
but can make a difference for individual clusters (see Section 
EI. We then stack the WMAP7 data by averaging ATj{9) 

over sources, j. 

We have followed ISawangwit fc Shanka (|2010al ) in using 
jack-knife errors, for both the radio and SZ sources, based 
on 6 equal area sub-fields defined by lines of constant galac- 
tic longitude and split by the galactic equator. We have also 



5 PLANCK RADIO SOURCE RESULTS 

5.1 Flux Comparisons 

We now compare WMAP7 sources at Q, V, W to their 
counterparts in the Planck ERS at 100, 70 and 44 GHz. 
We also compare the Planck fiuxes in the lOOGHz band to 
the gro und-based ATCA and IRAM s ource fiuxes previously 
used bvlSawangwi t fc Shanka (|2010al ). 

In Fig. [T] we compare the WMAP7 fiuxes to Planck and 
the ground-based ATCA and IRAM sources. We only con- 
sider the matches with separation less than 2' to avoid any 
possible systematic errors associated with sources that have 
poor astrometry. However, our results are independent of 
this cut up to 10'. At high fiuxes we see evidence for a sys- 
tematically lower WMAP flux, « 50% above 2Jy. This non- 
linearity is particularly prominent in the W band, the band 
with the greatest angular resolution. This Planck- WMAP 
non-linearity is statistically significant in W and V and ro- 
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1 10 

PLANCK S„,KH.[Jyl 



Figure 2. A comparison of the lOOGHz Planck fluxes and the 
ground-based sources. 



bust to high flux outliers. We find best fit logarithmic slopes 
of [0.70 ±0.04, 0.83 ±0.04, 0.95 ±0.05] for [W,V,Q]. Thus the 
one-to-one relation is strongly rejected in both W and V, al- 
though not in Q. 

Clearly, Planck and WMAP fluxes for sources were 
measured at different times. Since at least ~ 30% of the 
WMAP5 rad i o sou rces exhibit some level of variability 
IWright et al.l 12003 ) , we expect and observe much larger 
scatter than accounted for by the estimated flux uncertainty. 
We note that since the brighter WMAP sources are too faint, 
this is in the opposite sense expected if variability was bias- 
ing faint Planck sources into the WMAP catalogue when in 
a bright phase. 

The non-linearity between Planck and WMAP is in 
contrast to direct comparisons between Planck and ground- 
based data. These instead show good agreement, as shown in 
Fig. [2] for the Planck lOOGHz radio point sources. The best 
fit logarithmic slope of [0.97 ±0.05] is statistically consistent 
with the one-to-one relation. This is particularly significant 
because this is the frequency band with the greatest Planck- 
WMAP non-linearity. 

Given this agreement between Planck and the ground- 
based observations, we interpret the fiux non-linearity as 
being due to WMAP. In Section[7]we shall relate this prob- 
lem to the wider t han expected beam profiles observed by 
ISawangwit fc Shanks, (,2010a '). 



5.2 Planck Profiles 



We now repeat the stacking analysis of lSawangwit fc Shankd 
l|2010al ) centring on the Planck radio point sources. Planck 
sources are selected at multiple wavebands which may 
be advantageous in avoiding spurious sources etc. Figs. [3] 
(a),(b),(c) are shown for completeness because these raw 
temperature plots demonstrate the main uncertainty in this 
analysis which is the accuracy of the background subtrac- 
tion. We note that there is some difference between the 
global background (dotted line) and the background local to 
the source samples but generally this effect appears smaller 
in the WMAP7 data (eg at W) than it was in the WMAP 5 
datasets used bv lSawangwit fc Shankd l|2010al ). 



Figs. [3] (d),(e),(f) show the same profiles now back- 
ground subtracted and scaled to unity at the origin to 
produce b (9). We have used the 'photometric' subtrac- 
tion to produce the radio point source temperature pro- 
file, ATradio{0). For the WMAP7 dataset there is very lit- 
tle difference in the profiles resulting from global or lo- 
cal/photometric background subtractions. 

These b (6) are now compared to the WMAP Jupiter 
beam and the best fit bea m to the bright WMAP ra - 
dio source profiles found by ISawangwit fc Shanka (|2010a|) 
(dashed orange line in their Fig. 2). There is again evidence 
that the Planck selected radio sources suggest a wider beam 
than the Jupiter beam, particularly in the W band, although 
the Planck sources l ie slight ly below the profile fits from 
ISawangwit fc Shanka l|2010ah . However, statistically, both 
results are in good agreement. 

The normalisation of b {9) to unity at small scales forms 
a further uncertainty in these beam comparisons. In Figs. [3] 
(g),(h),(i) we have applied the formalism of Section 3.1 and 
attempted to make absolute normalisations of the various 
model profiles, using the Planck ERCSC listed fiuxes. Here 
if the radio sources followed the Jupiter profile, for example, 
we should see the same peak temperature for the stacked 
model profile and the stacked data. We see that the Planck 
peak temperatures, particularly in the W band, tend to lie 
between the Jupiter profile and the previous WMAP bright 
radio source fits. These results suggest that the previous 
radio source fit may be too wide at 6 > 30' where it is 
essentially an extrapolation, unconstrained by the data, and 
this will affect the accuracy of its absolute normalisation i.e. 
there is a large error in fibGam- Otherwise, the conclusion is 
similar to that from Figs. [3] (d),(e),(f) in that the Planck 
data is suggesting that the Jupiter beam is a poor fit to the 
radio source profiles particularly at W. 



6 WMAP AND NVSS RADIO SOURCES 



6.1 WMAP7 radio sources 



We next repeat the analysis of lSawangwit fc Shan ks '2010a|) 
for the WMAP7 data where the S/N is slightly higher than 
in the WMAP5 dataset originally used by the se authors. 
We co ntinue to use the WMAP7 source list of iGold et al.l 
(|201ll '). The results are shown in Fig. [l] We see that 
the results reasonably agree with the power-laws fitted by 



ISawangwit fc Shanks 



(|2010a ). At large scales, the profiles 
may return to zero more uniformly than previously, making 
the differences between the photometric and global profile 
estimates more marginal than previously found. 



6.2 WMAP7 'CMB-free' radio sources 

We take the opportunity to make measurements of the beam 
profile based on other radio source samples. One suggestion 
has been that the wider than expecte d radio source profiles 
might be caused by Eddington bias. ISawangwit fc Shanka 
(2010aj) have already argued against this on the grounds 
that effectively we selected on the GB6 and PMN ground- 
based data, and this meant that the profiles might not be 
widened due to mis-identifying upward CMB fluctuations 
in the WMAP maps as CMB sources. We now check this 
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Figure 3. (a),(b),(c): The raw stacked WMAP7 [Q,V,W] temperature profiles for the Planck [100,70,44] GHz band sources with the 
global mean temperature of the map plotted as the dashed line. (d),(e),(f): The photometrically subtracted, stacked and re- normalised 
WMAP7 [Q,V,W] b^ie) profiles for the Planck [100,70,441 GHz band sources. Also shown are the b''{9) for the Jupiter beam (blue, solid) 
and the radio source fit (red, dashed) of Sawangwit & Shanks (2010a). (g),(h),(i): The photometrically subtracted stacked WMAP7 
[Q,V,W] AT{6) profiles for the Planck sources, now absolutely normalis ed via the Planck flux. Also shown are the AT{6) for the Jupiter 
beam (blue, solid) and the radio source flt (red, dashed) of Sawangwit fc Shanka (|2010a|). 



result further by using the CMB-free WMAP catalogue of 
iGold et al.l (|201ll ). 



In the 'CMB-free' method, WMAP sources are selected 
using the Q,V, W bands simultaneously to form an internal 
linear combination map (ILC) with weights chosen to can- 
cel out the CMB anisotropy signal. Again any Eddington 
bias due to CMB fluctuations should be reduced in the case 
of this new point source catalogue. We therefore repeated 
our stacking ana lysis with the 417 QVW sources from the 
iGold et al.|[201ll WMAP7 catalogue (see Fig. [5]). Overall we 
again see wider-than-expected profiles, broadly consistent 
with the results in Figs. [3l [4] and the previous WMAP5 re- 
sults in Fig. 2 of lSawangwit fc Shanks! l|2010al ). 



6.3 NVSS radio sources 

Point source catalogues made at significantly lower frequen- 
cies than the WMAP bands are also unlikely to be affected 
by Eddington bias, if identification is done independently 
of the WMAP7 point source catalogue. For example, point 
sources selected at 1.4 GHz will have Rayleigh- Jeans tem- 
perature « 4500 X higher than a source with similar flux den- 
sity selected at W-band (« 94 GHz), i.e. Trj oc ^boam^"^' 
whereas the rms Rayleigh- Jeans temperature due to the 
CMB fluctuations stays roughly the sa me between the two 
frequency bands, (IBennett et al.l 120031 ). Therefore, we now 
stack WMAP7 temperature data centred around the posi- 
tions of the 1156 51.4 > 1 Jy NVSS point sources. Fig. \6\ 
shows the resulting Ql, VI and Wl profiles. We see that 
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Figure 4. (a),(b),(c): The raw stacked WMAP7 [Q,V,W] temperature profiles for the WMAP7 sources of lGold et all l|201ll ') with the 
global mean temperature of the map plotted as the dashed li ne. (d),(e),( f '): Th e photometrically subtracted, stacked and re-normalised 
WMAP7 [Q,V,W] b^iO) profiles for the WM AP7 sources of iGold et al.! feoilh . Also shown are the b=(e) for the Jupiter beam (blue, 
solid) and the radio source fit (red, dashed) of lSawanew^it &: Shanka l l2010al V 



1.00 



■h 0.10 



0.01 




Radio beam 
Jupiter beam b 



M WMAP7-cmbfree Q(362) 



Q(a) - 




V(b) 




W(c) 



1 10 

9 (arcmin) 



100 



10 
I (arcmin) 



100 



10 100 

e (arcmin) 



Figure 5. (a),(b),(c): The photometrica lly subtracted, stacked and re-normalised WMAP7 [Q,V,W] b''{8) profiles for the CMB-free 
WMAP catalogues of Gold ct al. (2011). Also shown are the b'^{9) for the Jupiter beam (blue, solid) and the radio source fit (red, 
dashed) of iSawangwit &: Shanks (2010a; ). 



they are consistent with thos e measured using WMAP5 to - 
tal/bright sources in Fig. 2 of lSawangwit fc Shanks! (|2010al ). 
However, the profiles do not appear as wide as the WMAP5 
faintest subsample given that the average flux of the NVSS 
sample at WMAP bands is « 3x lower. 



Many of the NVSS sources are resolved into multiple 
components (jBlake fc Walj|2002l ). However, this is unlikely 
to cause the widening of the beam beyond S > 6'. Here, 
as a precautionary measure, we shall test the beam pro- 
file measured using the NVSS by excluding any source that 
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Figure 6. (a),(b),(c):The photometrically subtracted stacked WMAP7 fQ,V, W] b^jO) profiles for NVSS sc) urces. Also shown are the 
6^(0) for the Jupiter beam (blue, solid) and the radio source fit (red, dashed) of lSawangwit &: Shankj l|2010a|). 



has neighbouring source (s) within 1°. This extra condition 
reduces the number of Si. 4 > 1.0 Jy sources outside the 
WMAP7 'point source catalogue' mask to 972. The result- 
ing beam profiles are also shown in the Fig.[Sl We see that 
the beam profiles are in good agreement with the previous 
results. 



7 RELATING WMAP FLUX NON-LINEARITY 
AND RADIO SOURCE PROFILES 

Our first attempts to explain the wider than expected ra- 
dio source profiles focused on the possibility that there 
was a timing offset between the WM AP antenna poin ting 
and temperature data, as proposed bv lLiu fc Lil (|201l|). As 
well as causing effects at large scale due to a wrongly sub- 
tracted dipole, t his scan pattern offset would cau se a wider 
beam profile (seelMoss. Scott, fc Sigurdsonll2011^ ■ However, 



NVSS and Planck profiles previously discussed. We see that 
the non-linear Jupiter beam appears to explain the ra- 
dio source profiles very well over the full range of angles. 
This suggests there may be a direct connection between 
the fiux non-linearity seen in Figs, [l] (a),(b),(c) and the 
wide radio source beam profiles. We also note that this 
'non-linear beam' may even be starting to approach the 
so-called 'diy' beam reverse engineered to fit a particular 
cosmological model with a fi rst acoustic peak at I — 330 
dSawangwit fc Shanks|[2010bl ). 

Taking the WMAP estimate of the Jupiter beam as 
being accurate when at ~ 1200Jy, which is significantly 
brighter than even the brightest « lOJy radio source we 
have used, is clearly an issue for the approach we have taken. 
Nevertheless, it is striking how well the average radio source 
profiles are reproduced when this assumption is made. 



ISawangwit et all ( in prepf) finds that while such an effect 8 SZ RESULTS 



could help match the radio sources at 10' scales, it could 
not explain their profile at 30' unless the offset was signifi- 
cantly greater than 25. 6ms in W band (~ half a p ixel) which 
may be unreasonable (jSawangwit et alllin prepf) . Although 
this suggests that such an offset did not occur at the map- 
making stage, there is still the po ssibility that it could have 
occured at the calibration stage IjLiu fc Lill201ll : iRoukemal 
I2OIOI ). 

We now look instead for an explanation for the ra- 
dio source beam profiles based on the non-linear WMAP 
radio-source fluxes and the wide beam proflles. We therefore 
take the power-law fits to the relations in Figs. [T](a),(b,(c) 
and apply them to the symmetrised Jupiter beam b^ (0) 
profiles which we assume is the true optical profile of the 
WMAP beam. The rationale for doing this is explained in 
Appendix [Bl where we assume that it is due to a non- 
linearity in the detector gain, g — fc(Ar)~°. So, taking 
the W band as an example, we applied the transformation 
b^'ie) = b^^.teriOY'"'^ = &lp«ter(e)°'', based On Fig.IHa), 
to the W beam profile as measured from Jupiter to form the 
effective beam profile, b (6), for the W band. 

The W band result is shown in Fig. [7| where it is com- 
pared to the error weighted mean of the empirical WMAP, 



8.1 WMAP-Planck SZ comparison 

We next compare the stacked WMAP7 temperature pro- 
files for 159 clusters listed in the Planck ESZ catalogue. We 
are using the 'photometric' approach to background sub- 
traction, with an annulus from 60' — 120' being used in W 
(and scaled according to beamwidth in Q and V). The self- 
similar models supplied by Planck as fits to the Planck data 
have been convolved with the WMAP7 Jupiter beam in the 
first instance. We see in Fig. |8] that these models are poor 
fits to the WMAP data. In [W,V,Q] they he a factor of 
Ki [2.2, 2.4, 2.8] below the models at ^ 10', with the discrep- 
ancy getting worse as we move to the lowest frequency Q 
and V bands. We also show similarly convolved isothermal 
models with parameters as listed in Section Although the 
self-similar fits are lower than these at 60' scales and thu s 
less discrepant with the WMAP data (|Koniatsu et al.ll201ll ). 
the discrepancy of both the isothermal and self-similar mod- 
els with the WMAP data remain poor at smaller scales. Of 
course, the main model here for inter-comparing WMAP 
and Planck is the self-similar model because that is the 
model fitted by the Planck team to their data. 

We first tried to explain the SZ discrepancy via the 
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Figure 7. The weighted mean of the stacked radio source beams in the Q, V and W bands from WMAP, NVSS and Planck catalogues 
compared to the Jupiter beam (blue), and the 'non-hnear Jupiter beam' (red) produced from applying the n on-linear WMAP-Planck 
flux relation in Fig. [Jla) to the Jupiter beam. Also shown is the 'DIY' b eam of lSawanewit &: Shankj ll2010bl ) reverse engineered from 
the raw WMAP Ci and the low Hq, ^baryon = 1 model of lShankj lIlQSSl V 



beam. In Figs. [8] we show the Planck s elf-similar models con- 
volved witii the power-law beams from lSawanerwit fc Shanks! 
l|2010ah . We find that in the case of the W band where 
the radio source profiles are most different from the Jupiter 
beam, there is some improvement in the agreement with the 
WMAP data. But even here the deficit is reduced only by « 
20% in the centre. Indeed, at larger scales, the wider beams 
slightly worsen the agreement with WMAP. In the Q and 
V bands where the radio source profiles are closer to the 
Jupiter beam, the wider beams give virtually no improve- 
ment in the agreement with the WMAP data. 

It is also the case that if we apply the non-linear flux- 
flux relation to the SZ proflles as we previously did to the 
radio profiles, we again see only a limited improvement in the 
agreement with the data. Again the WMAP SZ discrepancy 
is largest in Q where the radio source flux-fiux non-linearity 
is smallest, arguing that a different explanation must apply 
to the WMAP SZ results. 



8.2 Coma 

We have also looked at the Planck model fits for the Coma 
cluster and compared them to WMAP (see Fig. |5|. Part 
of the motivation here is that any beam problem may be 
expected to show up less for a nearby and therefore very 
extended {Osr^oo ~ 200') cluster like Coma, than in the case 
of more distant clusters. But here we see that the Planck 
model fit for Coma is now substantially overestimated by 
the WMAP data. This again argues that the explanation 
for the Planck- WMAP SZ discrepancy lies away from the 

beam profile. 

Previously iLieu Mittaz. fc Zhari3 (|2006l ') and 
iBielbv fc Shanka ( 20071 ) used an isothermal model to 
fit the Coma cluster. They used parameters estimated from 
the X-ray observations of Coma and found that the WMAP 
data could be fitted with these parameters. We show that 
this is still the case for the WMAP7 data. The reason for 



the discrepa ncy with the self - simila r fits to the Planck data 
is unknown. iKomatsu et al.l (|201ll ) has discussed whether 
Coma may sit on a large downwards CMB fluctuation but 
clearly this cannot explain the WMAP-Planck discrepancy. 



8.3 Explanations for the WMAP-Planck SZ 
discrepancies 

First, we discuss whether the WMAP SZ discrepancy may 
instead be caused by degeneracies in the self similar model 
fits derived from the ESZ. One issue is that these fits are 
parameterised by the integrated quantity, y(5-R5oo), yet the 

WMAP discrepancy is most significant at 6 < 10' where 
the small solid angle reduces its contribution to y(5i?5oo). 
However, the Planck fit is made by scaling the self-similar 
profile vi a an empirically measured central Compton param- 
eter, yp (JMehn. Bartlett. fc DelabrouiUd [2OO6I : iMelin et all 
I2OIII ). Therefore j/o is directly observed by Planck, suggest- 
ing that the small angle discrepancy may not be explained 
by an ambiguity in the ESZ SZ template. 

Bielby & Shanks previously discussed the possibility 
that contamination of SZ clusters by radio sources could 
have reduced the WMAP SZ signal. Because Planck only 
use ^ 100 GHz frequencies to identify SZ clusters, it could 
still be possible that the much smaller WMAP decrements 
seen at lower frequencies, particularly at Q(41 GHz), could 
be explained by radio source contamination only seen in the 

WMAP bands. However, the arguments against contamina- 
tion previously applied by Bielby & Shanks to the W band 
almost equally apply to the V and Q bands; the contami- 
nating sources fou nd in low frequenc y (eg 30GHz) surveys 
by authors such as lCoble et al.l ((20071) seem to have too low 
fiuxes to explain the WMAP SZ discrepancy. In addition, 
when we exclude « 10 ESZ clusters that have a contami- 
nating 44, 70 or 100 GHz source within 1 degree, this was 
found to make little difference to our stacked SZ results in 
any of the WMAP bands. 
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Figure 8. (a),(b),(c): The stacked SZ decrements of the Planck sources for the [W,V,Q] bands wi th the self-similar and isother mal /3 
models convolved with a WMAP Jupiter beam and a beam fitted to the radio source profiles by ISawangwit fc Shanksl l|2010a|). The 
stacked WMAP SZ decrements are always significantly weaker than any of these Planck-hased models. 



We note that iBielbv fc Shanksl (|2007l ') already found 
a similar WMAP disc repancy with the SZ results of 
iBonamente et al.l (|2006l ) whose results are now weighted 
heavily to small scales because of the high resolution of 
these interferometric observations. We have checked a stack 
of 12 Planck SZ clusters that also belong to the Bonamente 
et al catalogue and find good agreement in their SZ decre- 
ments. These results argue that the WMAP SZ discrepancy 
aX 6 < 10' is replicated in comparisons with other experi- 
ments covering a wide range of scales, motivating us to look 
for an explanation in a systematic WMAP instrumentation 
error. 

The large discrepancy seen at Q strongly argues that 



any systematics may be unrelated to the non-linear fluxes 
and wider beams measured for the radio sources. Since any 
optical beam profile effect must apply to both positive and 
negative AT equally, a more likely instrumental explanation 
could centre on different WMAP non-linearities at positive 
and negative AT (see Section [9J. 



9 DISCUSSION 

T he main criticism that was m ade of the previous results 
of ISawangwit fc ShanQ l|2010ah was that the wide WMAP 
radio source profiles may be caused by Eddington bias 
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Figure 9. (a),(b),(c): The Coma SZ decrements for the [W,V,Q] 
bands with the self-similar and isothermal /3 models convolved 
with a W MAP Jupiter beam and a beam fitted to the radio source 
profiles by ISawangwit fc Shanka 1120103^ . The error within each 
annulus for this individual cluster is simply the standard devia- 
tion of all clusters and is therefore only indicative. The observed 
WMAP decrements are now stronger than the Planck models and 
we also show an isotherma l model based on the X-ray parameters 
used bv iLieu. Mittaz. fc Zhang (2006); Biclby & Shanks (2007), 
convolved with the radio profile beam and which fits the WMAP 
data well. 



l|Eddingtonll 19131 ). Essentially, low S/N sources detected in 
WMAP data itself may be contaminated by upwards CMB 
fluctuations and not balanced by downwards fluctuations. 
This could explain the wider than expected profiles, partic- 
ularly at faint fiuxes. 

There may be some evidence for Eddington bias in 
the fain test W band source sub-sam ple that was initially 
used bv ISawangwit fc Shanks! (|2010al ). and for this reason 
the faintest sources were not used in their fits of the beam 
profile. Nevertheless, Eddington bias cannot affect sources 
detected at frequencies where CMB fiuctuations are sub- 
dominant. We have found that NVSS selected sources at 
low frequency also show a wider beam at W. Also we note 
that the Planck sources show the wider beam independent of 
whether the CMBSUBTRACT flag applies. WMAP sources 
selected from a 'CMB-free' map also show the same wider 
than expected beam^ 

Furthermore, ISawangwit et all l|in preo) ran Monte 



Carlo re-simulations of the source detection, producing artif- 
ical source catalogues extracted from simulated CMB maps. 
Here, after applying the same cross-correlation technique as 
for the data, the WMAP beam was recovered as input, again 
arguing that these sources are little alfected by Eddington 
bias. These authors also argued that an offset in the WMAP 
scan pattern could contribute to widening the beam profile 
bu t only if the off sets were larger than previously dscussed 
bv lLiu fc Lil (|201ll ') and then the fit at the largest, 30', scales 
could only be improved at the expense of the small-scale fit. 

Importantly, the Planck data also confirms the non- 
linearity of WMAP fiuxes, particularly in the W band. De- 
creasingly non-linear effects are also seen at Q and V. We 
have fitted this non-linearity and then applied the fit to the 
Jupiter beam profile. This procedure produced a surprisingly 
good fit to the radio-source profiles. 

The interpretation of the comparison of Planck- WMAP 
SZ decrements is less clear. The possibility that this discrep- 
ancy is also due to a wider than expected WMAP beam 
profile is made less likely by the fact that the Q and V band 
discrepancies are at least as bad as at W, but the radio 
source profiles in these bands are much less discrepant than 
the profile at W. The WMAP SZ decrements are not only 
discrepant with the Planck decrements but a lso with the 
the OVRO SZ decrements of iBonamente et all {2006). This 
suggests that the discrepancy may not be explained away by 
any degeneracies in the self-similar models that parameterise 
the Planck SZ results in the ESZ. 

If the SZ discrepancy with Planck is real, then it is 
possible that there is a different gain dependence at negative 
AT than positive AT. The radio sources can only check the 
WMAP calibration at positive AT. Indeed, it is possible 
that the gain is linear although with a different gain constant 
given that dividing the Planck models by a factor of « 2 
appears to bring them into agreement with the WMAP SZ 
decrements. 

It should be possible to test the calibration at both pos- 
itive and negative AT by comparing the goodness of fit of 
a dipole model with a linear and non-linear gain with both 
the raw and calibrate d counts and this work is in progress 
(jMalik et al.l lin prepf ). A preference for a non- linear gain 
from the raw data would indicate non-linear detectors. Oth- 
erwise if the raw data showed a linear gain then it would 
suggest that any non-linearity had arisen during the com- 
bined calibration/map-making process. 

We finally show in Fig. [TO] the effect the 'non-linear 
beam', created via AT (e) a AT{ef-'^, (see Fig. [Tt) has on 
the WMAP W1/W2 Ci. We assume here, for the moment, 
that this non-linear beam applies to both positive and neg- 
ative AT. Then debeaming the raw W1 /W2 Ci from Pol- 
Spice (Szapudi. Prunet. fc Colombill2001^ via equns (1,2) of 
ISawangwit fc Shanks! ( 2010al ). we see not only that the de- 
beamed power spectrum (blue line in Fig. llOp is significantly 
higher in amplitude, but that the first peak has moved from 
I = 220 to I — 275 compared to the model debeamed us- 
ing the Jupiter beam (cya n line). This is a similar effect as 
was reverse engineered bv !Sawangwit fc Sh anks (2010b) to 
fit the low Ho, fib = 1 cosmological model of Shanks (1985|) 
where the peak occurs at I — 330. Here the black line rep- 
resents the model and the red line the C; resulting from 
debeaming with the 'diy' beam of Fig. [T] The fact that the 
WMAP Ci debeamed with the 'non-linear' Jupiter beam has 
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Figure 10. The WMAP W1/W2 raw Q is shown as the green 
hne. The usual C; result from debeaming with the Jupiter beam is 
shown as the cyan line. The result produced by debeaming with 
the 'non-linear Jupiter beam' from Fig. [7] is sh o wn as the blue 
line. The low Ho, ^baryon = 1 model of Shankj l ll98a) is shown 
as the black line, together with the C; revers e engineered from this 
model via the 'diy' beam, also from Fig. [7] llSawangwit &: Shankj 
[2010b|) (red line). Also shown are the I > 600 results from the 
QUAD experiment l lBrown et al.ll2009h . 



a first peak now approaching the position of this very sim- 
ple model emphasises the importance of de-beaming to the 
cosmological interpretation of the WMAP power spectrum. 



10 CONCLUSIONS 

We have investigated the beam profile of WMAP by com- 
paring beam profiles from radio sources with the Jupiter 
beam profile. We have compared sources from Planck, NVSS 
and WMAP CMB-free catalogues. We find that in all cases 
the radio sources show wider profiles than the Jupiter beam 
with little indication of Eddington bias or dependence on the 
method of normalisation. We have also found non-linearity 
in WMAP fiuxes compared to Planck and ground-based 
fiuxes. Moreover, we have now related the fiux non-linearity 
to a gain non-linearity that can also explain the wide WMAP 
beam, particularly at W. 

We have also compared stacked WMAP SZ decrements 
with those measured by Planck and by ground-based obser- 
vations. Again, we find discrepancies with WMAP showing 
much lower decrements than observed in the Planck or in 
the ground-based data. The discrepancy is as large at Q as 
it is at W, arguing against an explanation involving the wide 
radio source profiles. We have suggested that there may be 
a different non-linear gain or at least a different linear gain 
that applies at negative AT. 

We have shown that transforming the Jupiter beam us- 
ing the non-linear radio source fiux relation results in an 
significant change to the amplitude and position of even the 
first acoustic peak in the CMB power spectrum. At the least, 
a wider beam would imply a much larger uncertainty in the 
normalisation and hence the estimate of erg from WMAP. 
Any uncertainty in the position of the first acoustic peak im- 



mediately translates into increased uncertainty in the CMB 
evidence for the CDM cosmological models. Clearly it is 
important to better understand the calibration and beam 
profile of WMAP, particularly in the W band which has 
the highest spatial resolution. Further tests using the dipole 
are currently being made on the raw and calibrated WMAP 
TOD to see if there is any suggestion that a non-linear gain 
solution may bo preferred to a linear g ain solution for the 
WMAP detectors (jMalik et alllin pretJ) . 
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APPENDIX A: SZ SELF SIMILAR MODEL 

In the self-similar SZ model as employed in the Planck ESZ, 
the fundamental parameters of a cluster are P500, M 500 and 
-Rsoo- Using the terminology of I Arnaud et ahl l|20ld ). 



M500 

^500 



47r 3 

-5- ^500 500pcrit 



Daiz) 



S5R500 



(Al) 
(A2) 



A I500 parameter corresponding to these is also defined, 



^50 



o-t 



(A3) 



which can be used as a characteristic SZ parameter instead 
of -Psoo- In equn. (|A3|I the units of Y500 are Mpc^, but are 
easily convertible to the arcmin^ units used in ESZ and 
throughout this papeQ. This Y500 is a distinct quantity from 
Y{R5oo) as found by evaluating equn. (|9}, In fact both P500 
and Fsoo ultimately cancel out of our calculation of the in- 
tegral Ycyi{< R), from the Planck provided parameters, the 
fully integrated Y5H500 and Osrsoo ■ 

However, the i ntroduction of yspp is well motivated be- 
cause, as shown by lArnaud et al.l l|2010r) . it allows a scale- 
free description of Ysph and Ycyi in terms of a; = R/R500 as 
follows. 



Y^phix) = Y^ooHx), (A4) 

Y^y,{x) = Ysph{5R5oo) - Y5ooJ{x) (A5) 

where I{x) and J{x) are the spherical and cylindrical scaling 
functions, 



I{x) = / 3V{u)udu, (A6) 

Jo 

J{x)^ 2,V{u){u -x^f^'^udu. (A7) 

J X 

We therefore find that 



Y,yi{x) = Y^oo{I{^)-J{x)), (A8) 

and since Ycyi{i) — Ysph{5) ~ /(5)Y5oo we finally obtain 



v;, 



,(i?) = yc„K5i?50o)(l-^) (A9) 



APPENDIX B: FLUX NON-LINEARITY 

The WMAP team assumption is that the detector gain, g, is 
independent of temperature and then that Cram ~ gx AT+b 
where Craw is the raw differential count and b is the baseline. 
We now assume a non- linear gain g = fcAT"". Then if the 
raw detector counts, Craw, are c^am ~ g^^T, (subsuming the 
baseline, 6, into Craw) then the recovered AT assuming a 
linear gain is 



AT: 



gAT/g ^ AT = kAT^ 



Y[Mpc^ 



.{-^f(Da[Mpc]fY[arcmin'^] 



UATa{9) = AT" 6-^ (61) where AT°{e) is the central antenna 
temperature and b (9) is the true beam profile normalised 
to unity at the centre, then 
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(l-c) 



ATf-^h'^ie)^'-"^ 



ATy(e) 



ATa{e) ex ATa{e) 

Since 6^(0) = 1, AT^ = Ar°''"°', so 6*'(6») = 6^(6i)<'-"' 
i.e. if the gain is non-linear to the power, —a, then a point 
source profile will be non- linear to the power, (1 — q). 

Next we look at what the above non-linear gain rela- 
tion means for fluxes. The relation between central antenna 
temperature and flux, S^, is 



Ar« = |^^x ^ 



Qb dT g(i^) 



— Oi/L V 



Therefore S^, = AT^/V^, and S^, = AT^/T^ where l/T^ = 
^B^^f-giy)- By using the Jupiter beam to convert peak 
temperature to a flux 5^, WMAP have effectively assumed 
Q.B = ^B and therefore F = F and so. 



S./S. = AT2/AT2 = Ar« = S:VZ 

Hence S^ = S^^^^'f^". So a non-linear gain with a AT"" 
temperature dependence produces non-linearity in both the 
radio profiles and the fiux scale to the (1 — a) power. 



This paper has been typeset from a TJ5X/ I^T[5X file prepared 
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